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H : 1 Introduction 

c^ , 

Quantum affine symmetries [1] - [5] are intensively developed as rich mathematical structures 
linking various branches of mathematical physics, like topological quantum field theories, inte- 
grable lattice models of statistical physics, rational conformal theories. They are implemented 
with the ultimate goal to understand and explore the consequences of the symmetries for de- 
scription of physical systems. 

Our work concerns the application of quantum symmetries for the exact solvability of many 
particle lattice systems interacting with stochastic dynamics. 

The main idea of integrability of lattice systems (within the inverse scattering method [B]) 
is the existence of a family of commuting transfer matrices, depending on a spectral parameter. 
For quantum spin chains the transfer matrices give rise to infinitely many mutually commuting 
conservation laws. This is the abelian symmetry of the system. The infinitely many commuting 



conserved charges can be diagonalized simultaneously and their common eigenspace is finite 
dimensional in most cases. Thus the abelian symmetry reduces the degeneracies of the spectrum 
from infinite to finite which is the reason for integrability. In addition many systems possess 
nonabelian symmetries. They determine the R matrix operator, a solution of the Yang-Baxter 
equation, up to an overall scalar factor and are identified as the quantum bulk symmetries. In 
the presence of general boundaries the quantum symmetry, and the integrability of the model 
as well, are broken. However with suitably chosen boundary conditions [TJ [8] a remnant of 
the bulk symmetry may survive and the system possesses hidden boundary symmetries, which 
determine a K-matrix, a solution of a boundary Yang-Baxter equation and allow for the exact 
solvability. Such non local boundary symmetry charges were originally obtained for the sine 
Gordon model [9] and generalized to affine Toda field theories [10], and derived from spin 
chain point of view as commuting with the transfer matrix for a special choice of the boundary 
conditions [11] or analogously as the one boundary Temperley-Lieb algebra centralizer in the 
" nondiagonal" spin 1/2 representation [T2] . 



In this paper we consider an algebraic prescription to construct two operators, possessing a 
coideal property with respect to the quantum affine Uq{sl{2)). We show that these operators 
generate an Askey- Wilson algebra which thus turns to be a coideal subalgebra of the quantum 
affine Uq{sl{2)). (See also [13] and [TUdS] for previous discussions of these problems.) We argue 
that one can construct a i^-matrix in terms of the Askey- Wilson algebra generators, which 
satisfies a boundary Yang-Baxter equation (known as a refiection equation). As an example of 
an Askey- Wilson boundary symmetry we consider a model of nonequilibrium physics, the open 
asymmetric exclusion process with most general boundary conditions. This model is exactly 
solvable in the stationary state within the matrix product ansatz to stochastic dynamics and 
it can be shown that the boundary operators generate the Askey- Wilson algebra. The model 
is equivalent to the integrable spin 1/2 XXZ chain with most general boundary terms, whose 
bulk Hamiltonian (infinite chain) possesses the quantum affine symmetry Uq{sl{2)). 



2 The Quantum Affine Uq{sl{2)) 

In this section we recall the definition of the affine Uq{sl{2)) [2], [3l [16]. We fix a real number 
< g < 1 (in the general case q is complex) and we use the g-symbol in the form 

W = gl/2 _ g-1/2 = W,V2 (1) 

The quantum affine Uq{sl{2)) is defined as the associative algebra with a unit with generators 
Ef and q^\ z = 0, 1 in the Chevalley basis and defining relations 

q^^q-"' = q-^-q^^ = 1 (2) 
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together with the g-Serre relations 

(EffEf - miEtfEfEf + [3]EfEf{Etr - Ef{Eff = 0, z ^ j (4) 

where [3] = l + q + q~^. The element c = Hq + Hi is central and its value is the level of the affine 
Uq{sl2). The algebra is endowed with the structure of a Hopf algebra. Namely, the coproduct 
A, the counit e and the antipode S are defined as 

(5) 



(6) 

1 (7) 

Let ao denote the longest root and p is 1/2 the sum of positive roots. We consider the Uq{sl{2)) 
algebra with a scaling element d, defined by {d, ao) = 1 ( (, ) is the nondependence bilinear 
form on the Cartan subalgebra) and denote h = (ao, ao) + 2(p, ao). With a finite-dimensional 
representation vry of Uq{sl{2)) one associates the quantum i?-matrix i?^^(A) which acts in 
V ®V and satisfies the Yang-Baxter equation. One has also [d ® d,R\ = 0. The universal 
i?-matrix -R(A) is uniquely defined [17] by the first terms in its expansion in powers of the 
Chevalley generators of Uq{sl{2)) 

^ _ ^cm+d^c+Y.UH,m^ ( 1 ® 1 + (gi/2 - g-i/2) ^ qHi/^Et ® q'^'^'^Er + , 

\ i=0 

One can define an automorphism T\ 

TxH, = Hi, TxEf = \^^Et, z = 0,l (9) 

and put -R(A) = {T{\) ® id)R. For a fixed finite-dimensional representation (vr, V) of the 
quotient algebra Uq{sl{2)), obtained by setting c = 0, one has i?^^(A) = {■n®n)R{\). Following 
[T7] . one introduces the currents L^{\) G EndV ®Uq{sl{2)), given by L+ = {id®nv){R),L~ = 
{id eg) 7ry)(-R*), which are explicitly expressed in terms of the Chevalley generators of Uq{sl{2)). 
The operators L^{X) generate a Hopf algebra A{R) and their matrix coefficients generate an 
algebra Aq{R) C A{R). Let -^(A) denote the quantum current 

L(A)=L+(Ag)(L-(A)-i (10) 

with a finite Laurent series expansion 

L(A) = ^/y(r^)A-"-^ neZ (11) 



A theorem (by Reshetikhin and Semenov-Tian-Shansky [T7]) states that the element t(A) = 
trg{L{X)) hes in the center of the quotient algebra of A{R), obtained by setting c = —h. Hence 
from the explicit expressions of the currents L^ in terms of the Chevalley generators follows 
that t(A) is the generating function of the Casimir elements of the quotient algebra of Uq{sl{2)), 
obtained by setting c = —h. 

For our purposes we will need a slightly different realization of the algebra in terms of 
the Chevalley generators and following [13] we define a new basis in Uq{sl{2)) generated by 
Hi, Ql, Ql 

Let now u,u*,v,v* be some (complex) scalars. We denote 

U = uQ'o, U* = u*Ql, V = vQ'o, V* = v*Ql (13) 

Then we have 

qV^VU - q~^'^UV = (g^/2 _ g-i/2)^^ (14) 

q^'^UV* - q~^/^V*U = (gi/' - q^'^^)uv*q~' (15) 

gi/V*f/* - g-i/2^*y* _ (^1/2 _ g-l/2)^*^* (Ig) 

ql/2u*y _ q-l/2yiJ* = (^1/2 _ ^-l/2)^*^^-c (^7) 

The operators U,U*,V,V* satisfy the following g-Serre relations which are the direct conse- 
quence of (4) 

U^U* - [3]U^U*U + [3]UU*U^ + U*U^ = (18) 

U*^U - [3]U*^UU* + [3]U*UU*^ + UU*^ = (19) 

V^V* - [3]V^V*V + [3]VV*V^ + V*V^ = (20) 

V*^V - [3]V*WV* + [3]V*VV*^ + VV*^ = (21) 

We can now consider the linear combinations 

A = U + V, A* = U* + V* (22) 

It can be verified directly by using the g-commutation relations (14-17) between the operators 
U,V,U*,V* and the g-Serre relations (18-21) that the operators A, A* satisfy the following 
relations which are the defining relations of a tridiagonal algebra 

A^A* - [3]gA^A*A + [3]gAA*A^ + A*A^ = -uv{q - q-^f[A, A*] (23) 

A*'^A - [3]qA*'^AA* + [3]qA*AA*'^ + AA*^ = -u*v*{q - q~^f[A*, A] 

where [3]^ = g -|- g^^ + 1 and [X, Y] = XY — YX. In the following section we are going to 
consider a general realization of the Askey- Wilson (AW) algebra as a coideal subalgebra of the 
quantum affine Uq{sl{2)) and the related tridiagonal algebra. 



3 The Askey- Wilson and tridiagonal algebras 

We begin this section with some definitions [HI [19], [201 [2I]- Let V^ be a vector space with 
(in)finite positive dimension. A tridiagonal pair on V is an ordered pair A, A* where A : V -^ V 
and A* : V ^ V are hnear transformations (see Definition 1.1 on p. 2 in [21]) and is a Leonard 
pair if the conditions hold: 1. There exists a basis for V with respect to which the matrix 
representing A is diagonal, and the matrix representing A* is irreducible tridiagonal. 2. There 
exists a basis for V with respect to which the matrix representing A* is diagonal and the ma- 
trix representing A is irreducible tridiagonal. (A tridiagonal matrix is irreducible whenever all 
entries on the superdiagonal and subdiagonal are nonzero.) A tridiagonal pair A, A* on V^ is a 
Leonard pair if for each A, A* all eigenspaces are of dimension 1. 

Definition 1: For a Leonard pair A, A* on V there exists a sequence of scalars /?, 7, 7*, p, p* ,uj, 77, 77* 
such that 

-[3AA*A + A^A* + A*A^--l{A,A*}-pA* = -f*A^ + uA + r] (24) 

-(3A*AA* + A*^A + AA*^--f*{AA*}-p*A = -fA*^ + uA* + 7]* 

which is uniquely determined by the pair. The equations (24) are called the Askey- Wilson 

relations. 

Definition 2: For a tridiagonal pair A, A* on V there exists a sequence of scalars (3, 7, 7*, p, p*, 

such that 

[A,-(3AA*A + A^A* + A*A'^ --f{A,A*} - pA*] = (25) 

[A*,-f]A*AA* + A*^A + AA*'^--f*{A,A*}-p*A] = 

which is uniquely determined by the pair, if the dimension of V is at least four. It has been 
proved [22] for a Leonard pair that the two definitions are equivalent. 

Definition 3: A tridiagonal algebra, respectively the Askey- Wilson algebra, is an associative 
algebra with a unit generated by the pair A, A* subject to the relations (25), respectively (24). 

Affine transformations 

A^tA + c'. A* -^ t*A* + c* (26) 

where t,t*,c',c* are some scalars, act on the operators A, A* of the Askey- Wilson relations 
and can be used to bring a tridiagonal pair in a reduced form with 7 = 7* = 0. (The label 
c' is used to avoid confusion with the Uq{sl{2)) central element.) Examples are the g-Serre 
relations with /5 = g + q~^ and 7 = 7* = p = p* = and the Dolan-Grady relations [23] with 

/? = 2, 7 = 7* = 0, p = P, p* = k*^ 

[A,[A,[A,A*]]] = e[A,A*] (27) 

[A*, [A*, [A*, A]]] = k*^[A*,A] 

The algebra (24) was first considered by Zhedanov [24j who showed that the Askey- Wilson 
polynomials give raise to two infinite-dimensional matrices satisfying the AW relations. The 
tridiagonal relations have recently been discussed in a more general framework [TSl dSl [2S] and 



a particular case of a tridiagonal algebra with p = p* through a homomorphism to Uq{sl{2)) 
was proposed in [15]. Leonard pairs have been classified according to their dependence on the 
sequence of scalars p^ [T9] and a correspondence to the orthogonal polynomials in the Askey- 
Wilson scheme was given. In |2S] the AW algebra (24) with 7 = 7* = has been equivalently 
described as an algebra with two generators and with structure constants determined in terms of 
the elementary symmetric polynomials in four parameters a, b, c, d, abed 7^ g"^, m = 0, 1, 2...; g 7^ 
0,q''^l,k=l,2,.... 

We can now formulate the following statement which defines a homomorphism of the AW 
algebra to the quantized affine algebra Uq{sl{2)) (see [15] for the relation to tridiagonal algebra). 

Proposition I: Let u, u* , v, v*, k, k* be some scalars. The operators A, A* defined by 

A = uE+q-"''^^ + vE-q~"°/'^ + kq~"" (28) 

A* = u*E+q-^'/^ + v*E{q-"'/^ + k*q-"' 

(It is assumed that Ef in (28) are rescaled by ±(g^/^ — q"^^"^) according to (12)) and their 
g-commutator 

[A, A*]g = q^^^AA* - q-^^^A*A (29) 

form a closed linear algebra, the Askey- Wilson algebra 

[[A,A*]„A], = -pA*-LuA-7] (30) 

[A*,[A,A*],], = -p*A-iuA*-7]* 

where the (representation dependent) structure constants are given by 

- p = uv{q-q-^Y, -p* = u*v*{q-q-^f (31) 

_ ^ = _(^V2 _ ^-1/2)2 ^j^j^. ^ ll[uu*q^'^ + v*vq-^/^)) (32) 

_ ry = (g _ g-i)(gV2 _ ^-1/2) (_fc(^^*gi/2 + ^*,,g-i/2) _ l^uvk*) (33) 

_ ry* = (g _ g-i)(gV2 _ ^-1/2) (_r (nn*^^/^ ^ ^*,,^-i/2) _ f^u*v*k) (34) 

The element ly is the coefficient lv{n) to A^"~^, for n = 0, in the Laurent series (11) of the 
quantum current for any highest weight module over Uq{sl{2)) and is central (the quadratic 
Casimir element) in the quotient algebra of Ug{sl{2)), obtained by setting c = —h (with h = 1). 
On a highest weight irreducible representation V ly is a. scalar. 

We note that eq.(28) defines the homomorphism of the AW algebra to the affine algebra 
Uq{sl{2)). For k = u + v,k* = u* + v* one recovers the particular case (22). Special cases of 
this homomorphism is the representation considered by Terwilliger [TS] with p = p* = and 
the one by Baseilhac [T^ [T5] and Baseilhac and Koizumi [SB] with u = v* and v = u* in (28). 
Making use of the evaluation representation for the Uq{sl{2)) generators in (28) 

7i,{Et) = E^, n,{E^) = u^'E^, n,(q^^) = q"" , T,,{q"^^) = q~^ (35) 

where E^,H are the Uq{sl(2)) generators, we obtain the Granovskii and Zhedanov realization 



The algebraic relations (23) follow from (30) by taking the commutator with A and A* 
respectively. The homomorphism (28) defines the AW algebra with two generators as the linear 
covariance algebra for the Uq{sl{2)) with operator valued structure constants. The tridiagonal 
algebra (23), denoted TD, is obtained through the chain of homomorphisms TD —>■ AW -^ 

U,{sl{2)). 

Definition 4'- The AW algebra with two generators A, A* defined by the homomorphism (28) 
as a coideal subalgebra yields a deformation in two parameters p, p* of the g-Serre relations of 
level zero quantum affine Ug{sl{2)), such that it results in a shift of the central charge to a non 
zero value c = — 1. 

Taking the limit q —>■ 1 in the defining relations (23) one obtains a two parameter deforma- 
tion of the Serre relations of level zero affine s/(2), known as the Dolan-Grady relations. 

From the explicit realization of the operators A, A* it follows that they generate a linear co- 
variance algebra for the Uq{sl{2)) which has the property of a coideal subalgebra. Let Bq{sl{2)) 
denote the algebra generated by A, A* . 

Proposition II: The Askey- Wilson algebra defined by the homomorphism (28) is a coideal 
subalgebra of Uq{sl{2)). The proof is straightforward by using the comultiplication (5). One 
has 

A(A) = /(g)A + (A-A;/)®g-^" (36) 

A(A*) = I®A* + {A*- k*I) ® q-"^ 

where the expressions on the RHS of (36) obviously belong to Bg{sl{2)) ® Uq{sl{2)). 



4 Representations of the Askey- Wilson algebra 

The Askey- Wilson algebra is known to possess very important properties which allow to obtain 
its ladder representations. We briefiy sketch these properties (for details see [HI [2H [25]). 
Namely, there is a representation with basis /,. with respect to which the operator A is diagonal 

Afr = Xrfr (37) 

where the eigenvalues satisfy a quadratic equation 

Xl + \1-{q + q-')\r\s-p = 0, (38) 

(which yields two different eigenvalues A^+i and A^._i for a fixed eigenvalue A^) and the operator 
A* is tridiagonal 

A* fr = a^+l/r+l -|- Kfr + C^-l/r-l (39) 

Depending on the sign of p the spectrum of the diagonal operator is hyperbolic of the form sh 
or ch and exp if p = 0. The algebra possesses a duality property. Due to the duality property 



the dual basis exists with respect to which the operator A* is diagonal and the operator A is 
tri diagonal. We have 



A*f; = a;/; 

Afs = (^s+lfs+1 +bsfs + Cs_l/s-l 



(40) 



where A* satisfies the quadratic equation (38) with —p replaced by —p*. The overlap function 
of the two basis {s\r) = {fs\fr) can be expressed in terms of the Askey- Wilson polynomials. 
Let Pn = Pnix;a,b,c,d) denote the n-th Askey- Wilson polynomial [28] depending on four 
parameters a, b, c, d 



Pn{x]a,b,c,d) =4 $3 



q "',abcdq"- ^,ay,ay ^ 
ab, ac, ad 



\q;Q 



(41) 



with pQ = l,x = y + y ^ and < g < 1. Then there is a basic representation of the AW algebra 
[25] which is a representation of the tridiagonal algebra as well [19] 

[A, A^A* - (g + q-')AA*A + A*A^ + abcdq-\q - q'^fA*] = 

[A*, A*^A - (g + q-^)A*AA* + AA*^ + (g - q-^fA] = 
in the space of symmetric Laurent polynomials f[y] = fly"^] with a basis {po,pi, ...) as follows 

Af[y] = {y + y-')f[y], A*f[y] = Vf[y] (42) 

where V is the second order g-difference operator [28] having the Askey- Wilson polynomials pn 
as eigenfunctions. It is a linear transformation given by 



^f[y] 



abcdq ^)f[y] + 



ay)il -by){l -cy){l - dy) 



+ 



{I - y^){l - qy^) 
[a - y){b - y){c - y){d - y) 



if[qy]-f[y]) (43) 



(/[rt]-/M) 



(1 -y2)(g-y2) 

with T>{1) = 1 + abcdq~^. The eigenvalue equation for the joint eigenfunctions pn reads 

T^Pn = XnPn, K = ?"" + abcdq""-^ 



(44) 



and the operator A* is represented by an infinite-dimensional matrix diag(Ao, A*, A2, ...). The 
operator Apn = xpn is represented by a tridiagonal matrix. Let A denote the matrix whose 
matrix elements enter the three-term recurrence relation for the Askey- Wilson polynomials 



XPn = bnPn+l + dnPn + CnP 



n—lj 



P-i = 



A 



fao ci \ 

bo ai C2 

bi 02 

V ■ 7 



(45) 
(46) 
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The explicit form of the matrix elements of A reads 

an = a + a~^ - b^ - Cn (47) 

_ (1 - a6g")(l - acg")(l - arfg")(l - a6cc/g"-i) 

_ a(l - g")(l - 6cg"-i)(l - Mg"-^)(1 - crfg""^) 

The basis is orthogonal with the orthogonality condition for the Askey- Wilson polynomials [2H] 

Pn = a~"'{ab, ac, ad; q)nPn 

1 

/wix) 
- — ^ Pm{x] a, b, c, d\q)Pn{x] a, b, c, d\q)dx = hJmn (50) 

ZTTVl — X^ 
-1 

Where w{x) = ''^X::^^'!^^^/'^' ^^^ ^(^' /^) = Hr^otl - 2/ixg'= + /.^g^^], and 

, _ [abcdq'''^; q)n{abcdq'^''^; q)oo , , 

(g"+^, abq"^, acq"-, adq", bcq", bdq"-, cdq"^; q)oo 

As noted in the previous section a tridiagonal pair of operators A, A* is determined up to 
the affine transformation. One can appropriately rescale the operators to obtain the algebraic 
relations in the form 

[A, A^A* - pAA*A + A*A^ - (g - q-^fA*] = (52) 

[A*, A*'^A - f3A*AA* + AA*'^ - (g - g~^)M] = 

which results in the corresponding rescaling of the matrix elements. A shift of the operators 
has no other effect but shifting the diagonal elements of the representing matrices. 



5 Askey- Wilson algebra and reflection equation 

We consider models of statistical physics in which the spin variable is associated with the site 
i of a one-dimensional lattice. An example of a model with quantum affine symmetry is the 
spin 1/2 XXZ model with Hamiltonian defined on an infinite dimensional chain 

i 

where the Pauli matrices af,a^,a^ act on the i-th component of the infinite tensor product 
... ® Vi-i ® Vi ® Vi+i ® ..., with V = C^. This model is known to be integrable j29] within the 



representation theory of the affine quantized algebra Uq{sl{2)). Namely, given the Uq{sl{2)) R- 
matrix operator R{zi/ Z2) G EndcVz^®Vz^, where Vz is the two-dimensional Uq{sl{2)) evaluation 
module, satisfying the Yang-Baxter equation 

Rl2{Zl/z2)Rn{zi)R2z{z2) = R23{Z2)R13{ZI)R12{ZI/Z2) (54) 

then, the Hamiltonian is written as if = ^ifjj+i, where the two site Hamiltonian density is 
obtained as 

Hii+i = -1- P Rii+i\u=o (55) 

with P the permutation operator and Zi/z2 = e". The generators act on the quantum space by 
means of the infinite coproduct and the invariance with respect to the afiine Uq{sl{2)) manifests 
in the property 

[i7,A-(Gfe)] = (56) 

for any of the generators Gk of Uq{sl{2)). If we introduce for finite chain a boundary of a 
particular form, such as diagonal boundary terms, the symmetry is reduced to Uq{sl{2)) and 
the invariant Hamiltonian has the form [30] 

hTx, = -1/2 Y.{aM+, + a^al, + A,a,Vf+, + hia^, - af) + A,) (57) 

where 

A, = -^(g + g^'), h = ^{q-q^') (58) 

In the presence of a boundary in addition to the i?-matrix there is one more matrix K{z) which 
satisfies the boundary Yang-Baxter equation, also known as a reflection equation. 

R{z,/z2){K{zi) ® I)R{ziZ2){I ® K{z2)) - (/ ® K{z2))R{ziZ2){K{zi) ® I)R{zi/z2) = (59) 

Within the quantum inverse scattering method the K-matrix is related to the quantum current 
L = L~^{L~)~^ (10) where L"^ G EndV ® Uq{sl{2)). In section (2) the two generators of the 
Askey- Wilson algebra were constructed as linear covariant objects with the coproduct properties 
of two-sided coideals of the quantum affine symmetry Uq{sl{2)). It is suggestive to construct 
the if-matrix in terms of the AW algebra generators. 

Let R{z) be the symmetric trigonometric i?-matrix with deformation parameter g^/^ 

/gl/2^ _ g-l/2^-1 \ 

p. . z-z-^ qy^ - g-V2 q ,„„. 

^(^) = n .1/2 .-1/2 . ^-1 n (60) 



acting on the auxiliary tensor product space Vz-^ ® Vz^ which carry the fundamental represen- 
tations of the covariance algebra. Then one can construct an operator L{z) [1] in terms of the 

Uq{sl{2)) generators 

^qJz _ z-^q-Jz (gi/2 _ g-i/2) j_ 



(q^^z - 


- q-^l^z' 


-1 





\ 











qi/2 _ q^i/2 









qi/2 _ q-i/2 


z-z-^ 


\ 










q'/^z - q~^'^z-^l 



L{z) - (^^^1/2 _ ^-1/2) j^ ^^-J3 _ ^^iqJ, ) (61) 
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acting on the tensor product Vq ® Vq of the auxihary space Vq and the quantum space Vq 
which in the general case carry finite dimensional inequivalent Ug{sl{2)) representations. The 
L-operator satisfies 

R{z,/z2)L,{z,)L2iz2) = L2{z2)L^{zi)Riz,/z2) (62) 

where Li = L®/ and L2 = I®L. As it is known this relation together with the refiection equa- 
tion (59) constitute the basic algebraic relations of the inverse scattering method to integrable 
models. 

We are now going to construct a solution to eq.(59) in terms of the operators A, A* . 

Proposition III: Let A, A* generate the AW algebra, the linear covariance algebra for Uq{sl(2)). 
Then there exists a refiection matrix K{z) = K°^{z) + K'^{z), constructed in terms of the AW 
algebra generators, where the part K°^ has the form 



q^ 



'/^zA - q-^/^z-'^A* -^(g^/2 _ q~^/^)[A\ A],' 



^ (^) - ^-lVPr^.1/2 ^-l/2\l A 4*1 ^-1/2^-1 A , ..1/2^ Vp a* ^^^^ 



i-^(gV2 _ q-y^)[A, A*], -q-'/'z-'A + q'^'z^A 



and the part K'^{z) is 



-P 



q^''^zr]* — q ^''^z ^r] q ^'"^ztj — q^''^z ^rj* 



K" = ^ -1 ^ - 1- K'' = - — - - — (64) 



_ gi/2^2 ^ g-1/2^-2 ^ ^ gi/2^2 ^ ^-1/2^-2 _ _^ ^ 

^12- P ^1/2 + ^-1/2 ^'^' ^21- ^1/2 + ^-1/2 P ^'^ 

The matrix K{z) is a solution of the boundary Yang-Baxter equation (59) provided the oper- 
ators A, A* obey the tridiagonal algebraic relations of the AW algebra in the reduced general 
form (30) with all structure constants p, p*,u!,ri,ri* nonzero. We denote this solution K{z,p). 

The proof of this proposition is rather long but straightforward. It is directly verified using 
the explicit form of the R matrix (60) and the AW algebraic relations (30) that the boundary 
matrix K from (63, 64) solves the refiection equation (59). 

We emphasize on the factor -^ to A* and to in the K matrix. This factor is due to the 
fact that the solution of the boundary Yang-Baxter equation (59) in terms of the AW algebra 
generators requires p = p*. This is not a problem since given the AW algebra in the general 
form (30) we can relate it to an algebra with p = p* rescaling A* -^ -^A*. Alternatively we 



can rescale A -^ -^A to obtain an AW algebra with p = p*. This gives a second solution 
K{z,p*) of the refiection equation. Its K°P{z,p*) part has the form 



i/2^VZ^_ -1/2 -1^* _VKfql/2_q^l/2^lJ^* M\ 

y) l_p*-i^(gi/2_^-i/2)[^^^*]^ _q~i/2^^i^^^qi/2^^*j y > 

The matrix elements of the K'^{z,p*) part are obtained from (64) by the interchange p *-> p*. 
The solution K{z,p*) can be implemented to construct a solution K*{z) of the dual refiection 
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equation [3, [31] Namely, the matrix K*{z) = K^{z ^,p*) solves the dual reflection equation, 
(which is obtained from eq.(59) by changing 2:1 2 -^ Q~^^'^z^2 ^^^ ^ ~* K^)- 

Setting p = p* and 77 = 77* = in (63) and (64) we obtain the K matrix considered in 
[?] for such a very particular case of an AW algebra and for the spin 1/2 quantum space 
representation. We note that an AW algebra in the general form with a sequence of scalars 

— (g+g~^), 7, 7*, w, ?7, r]* cannot be reduced to such a particular algebra with structure constants 

— {q + q~^), p, p,0,0,uj,0,0. There exists an unique affine transformation [52] to only set 7 = 
7* = (and simultaneously either p = r]* = or r] = p* = 0). 



6 A model of nonequilibrium physics with boundary Askey- 
Wilson algebra 



Reaction-diffusion processes provide a playground to increase the utility of quantum groups 

As a physical example we consider the asymmetric simple exclusion process (ASEP), a model 

of nonequilibrium physics with rich behaviour and wide range of applicability [3i 



The asymmetric exclusion process is an exactly solvable model of a lattice diffusion system 
of particles interacting with a hard core exclusion, i.e. the lattice site can be either empty 
or occupied by a particle. As a stochastic process it is described in terms of a probability 
distribution P{si,t) of a stochastic variable Sj = 0, 1 at a site i = 1,2, ....L of a linear chain. 
A state on the lattice at a time t is determined by the occupation numbers Sj and a transition 
to another configuration s'^ during an infinitesimal time step dt is given by the probability 
T{s,s')dt. Due to probability conservation r(s, s) = — X^^/.^^ r(s', s). The rates T = F*;, 
i,j,k,l = 0,1 are assumed to be independent from the position in the bulk. For diffusion 
processes the transition rate matrix becomes simply F^^ = Qi^. At the boundaries, i.e. sites 
1 and L additional processes can take place with rates L^ and Rl {i,j = 0, 1). In the set of 
occupation numbers (si, S2, ..., sl) specifying a configuration of the system Sj = if a site i is 
empty, Sj = 1 if there is a particle at a site i. Particles hop to the left with probability goidt 
and to the right with probability giodt. The event of exchange happens if out of two adjacent 
sites one is a vacancy and the other is occupied by a particle. The symmetric simple exclusion 
process is known as the lattice gas model of particles hopping between nearest-neighbour sites 
with a constant rate qqi = qiq = g. The partially asymmetric simple exclusion process with 
hopping in a preferred direction is the driven diffusive lattice gas of particles moving under 
the action of an external field. The process is totally asymmetric if all jumps occur in one 
direction only, and partially asymmetric if there is a different non-zero probability of both left 
and right hopping. The number of particles in the bulk is conserved and this is the case of 
periodic boundary conditions. In the case of open systems, the lattice gas is coupled to external 
reservoirs of particles of fixed density. Phase transitions inducing boundary processes [3B] are 
the most interesting examples (see [39] for a review) when a particle is added with probability 
adt and/or removed with probability '-/dt at the left end of the chain, and it is removed with 
probability j3dt and/or added with probability 6dt at the right end of the chain. Without loss of 



12 



generality we can choose the right probabihty rate qiq = 1 and the left probability rate goi = q. 

The time evolution of the model is governed by the master equation for the probability 
distribution of the stochastic system 

^-^ = J2T{s,s')P{s\t) (66) 

s' 

It can be mapped to a Schroedinger equation in imaginary time for a quantum Hamiltonian 
with nearest-neighbour interaction in the bulk and single-site boundary terms 

^ = -HP{t) (67) 

wheieH = 'YljHjj+i + H'^^^ + H^^\ The ground state of this, in general non-Hermitian, 
Hamiltonian corresponds to the stationary probability distribution of the stochastic dynamics. 
The mapping provides a connection to the integrable S't/q(2)-symmetric XXZ quantum spin 
chain with anisotropy A = ^^^ — -, q = ^^ ^ 1 and most general boundary terms. 

We consider the model within the matrix-product-state ansatz of stochastic dynamics [39l 
HO], which was inspired by the inverse scattering method to integrable models. The idea is 
that one associates with an occupation number Si at a position i a matrix D^- = Di ii a. site 
i = 1,2,...,L is occupied and Dg. = Dq ii a. site i is empty and the stationary probability 
distribution is expressed as a product of (or a trace over) matrices that form a representation 
of a quadratic algebra 

DiDq — qDoDi = xqDi — DqXi, xq + xi = (68) 

where < g < 1 and xcXi are representation dependent parameters. The totally asymmetric 
process corresponds to g = and the symmetric process - to g = 1. The quadratic algebra with 
no x-terms on the RHS (i.e. -Di-Do ^ qD^Di = 0) corresponds to a bulk process with reflecting 
boundaries. 

For an open system with boundary processes the normalized steady weight of a given con- 
figuration is expressed as a matrix element in an auxiliary vector space 

P(.......,.)^ <"'l^-^'-^--l"> . (69) 

with respect to the vectors \v > and < w\, defined by the boundary conditions 

{(3D,-6Do)\v) = xo\v) (70) 

{w\{aDo - -fDi) = {w\{-xi) 

The normalization factor to the stationary probability distribution is 

ZL = {w\{Do + Dif\v) (71) 
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In the following we set xi = —xq. 

The advantage of the matrix-product ansatz is that once the representation of the diffusion 
algebra and the boundary vectors \v > and < w\ are known, one can evaluate all the rele- 
vant physical quantities such as the mean density at a site i, (s,) = ^^'^ ""*" ^' — ^^^ "'^ ^' — —, 

the two-point correlation function (siSj) = ° — ^^^-^ — ^^ ^^^-2 — — and higher 

correlation functions. The current J through a bond between site i and site i -|- 1, J = 
HiDo+D^y-HD^Do^gDoD^)iDo+D,)^—^v) ^^^ ^ ^^^^ ^^^^^^ ^^^^ J = ^P^. The matrix Do + D, 

enters all the expressions and plays the role of a transfer matrix operator. 

The algebraic matrix state approach (MPA) is the equivalent formulation of recursion rela- 
tions derived for the asymmetric exclusion process (ASEP) in earlier works [1T],|12] which could 
not be readily generalized to other models. In most applications one uses infinite dimensional 
representations of the quadratic algebra. Finite dimensional representations [l3l HI] impose a 
constraint on the model parameters. They may be useful in relation to Bethe Ansatz on a ring 
|45j . The MPA was generalized to many-species models [391 HSl HZ] and to dynamical MPA 



The matrices Do, Di of the MPA generate an AW algebra with p = p* = [49]. We call this 
algebra the bulk Askey- Wilson algebra. For the particular case of only incoming (outgoing) par- 
ticle at left (right) end of the chain the boundary operators satisfy an isomorphic p = p* = AW 
algebra which can be solved by shifted g-deformed oscillators [50l ISTj [52] as they were applied 
for the ASEP with such particular boundary conditions [53l ^^ . In the general case of incom- 
ing and outgoing particles at both boundaries there are four operators /5-Di, —6Dq, — 7-Di, aDo 
and one needs an addition rule to form two linear independent boundary operators acting on 
the dual boundary vectors. From the quadratic algebra (68) two independent relations for the 
boundary operators follow 

pDiaDo - qaDopDi = xo{af3Di + f3aDo) (72) 

and 

-iDi5DQ - q6Do-fDi = xoi5-fDi + -fSDo) (73) 

To form two linearly independent operators B^ = j3Di—6Dq, B^ = —'jDi+aDQ for a solution to 
the boundary problem and in order to emphasize the equivalence of the ASEP to the integrable 
spin 1/2 XXZ one can use the Uq{su{2)) algebra. It is generated by three elements with the 
defining commutation relations 

[N, A^] = ±A± [A+, A.] = J,,~l,,, (74) 

and a central element 

The relations (72), (73) can be solved by choosing a representation of the boundary operators 
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[A+, 


,A. 


= 


g^- 


q~ 


-N 


gV2_ 


q' 


-1/2 


g^- 


1/2. 


-Q' 


-N+l/2 






{<! 


L/2_ 


-T 


-1/2)2 





in the form 



(3D,-6D, = (76) 

^0^ /V/2 . ^0^ A N/2 -(3q^'^ + 5 N (3-5 

aDo — '-fDi = 



VI - 9 VI - 9 1-9 1-9 

Separating the shift parts from the boundary operators and denoting the corresponding rest 
operator parts by A and A* we write the left and right boundary operators in the form 

/3D,-5Do = A + xo^ (77) 

1-9 

aZ?o - 7-Di = A* + xq- 

1-9 

Then the operators A and A* defined by 

A = pDi-6Do-xo^^ (78) 

1-9 

A* = aDo - -jDi - xq- 

1-9 

and their g-commutator 

[A, A*]q = q^l'^AA* - q~^'^A*A (79) 

satisfy the boundary Askey- Wilson algebra of the open ASEP 

[[A,A\,A], = -pA*-u;A-r, (80) 

[A\[AA*],], = -p*A-ujA*-7]* 

with structure constants given by 

p = xl(36q-\q'/' + q-'l')\ p* = xla^q'\q'/' + q-'/'f (81) 

-u; = xliP - 5)(7 - a) - x^Pj + a5){q''^ - q^^/^)Q (82) 

{[3 - 5){(3-t + a5) 



V 



q^q^^' + 9-/^)x^ (a,iP - 5)Q + (^ " ^ + f 7) 



gl/2 _ g-1/2 



3) 



Qri/2 _ g-1/2 



One can further use the afiine transformation properties of the AW algebra generators to 
obtain a representation of the boundary ASEP operators from the basic representation of the 
AW algebra. We summarize the results for the representation of the ASEP boundary operators 
(for details see [19] ): 
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There is a representation tt in a space with basis the AW polynomials p„ = Pn{x] a, b, c, d) 
(41) 

(po,Pi,P2,-)* (84) 

with respect to which the right boundary operator Di — -^Dq = Di + bdDo is diagonal. The 
representing matrix is diag(Ao, Ai, A2, ...) with the eigenvalues A„ given by 

A„ = -^ (6g~" + dg" + 1 + bd) (85) 

1 — q ' 

The left boundary operator Z^o — ^-Di = -Dq + acDi is tridiagonal and its representing matrix 
has the form 

^1/2 

Mq + acDx) = (bA* + 1 + ac) (86) 

1 — g ^ 

where the matrix A is given by (46). The dual representation it* has a basis 

iPo,Pi,P2,--) (87) 

with respect to which the left boundary operator 7r*(Do + acDi) is diagonal diag{XQ, X^, ...) 
with diagonal elements 

a; = ^ (ag-" + eg" + 1 + ac) (88) 

1 — g ^ 

The right boundary operator is represented by a tridiagonal matrix 

7r*(Di + MDo) = -^ (aA+l + bd) (89) 

1-g 

The ASEP boundary value problem is satisfied with the left and right boundary vectors 
chosen in the form 

{w\ = h',^'\p,, 0, 0, ....) \v) = h-^'\p^, 0, 0, ...)* (90) 

where h^ is a normalization from the orthogonality condition (50). With this choice the so- 
lutions to the boundary eigenvalue equations uniquely relate (in this representation) the four 
parameters of the Askey- Wilson polynomials with the boundary probability rates 

a = K^, 6 = K+, c= K*_, d = K_ (91) 

where 



^ ~i/^ - "^ - \^ - H) ) ^ y yh' - ^ - y'- - H) ) t ^i-^^ ,^^. 



-ip- 


-6- 


-(1- 


-q))± 


ViP- 


-5- 


-(1- 


-q)? 


+ 4:f3S 




-7- 


-(1- 


-q))± 


W 










-(«- 


V(«- 


-7- 


-(1- 


-<i)y 


' + 4a7 



^ 2a 



The expressions on the RHS of eq.(92) are the functions of the parameters which define the 
phase diagram of the ASEP. They have been used in previously known MPA applications 
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where have always been taken for granted. It is quite remarkable that here they follow from 
the properties of the Askey- Wilson algebra representations. 

It can be further shown that the transfer matrix Dq + Di and each of the boundary operators 
generate isomorphic AW algebras [19]. In the tridiagonal representation the transfer matrix 
Do + Di satisfies the three-term recurrence relation of the AW polynomials which was explored 
in [55] for the solution of the ASEP in the stationary state. The exact calculation of all the 
physical quantities, such as the current, correlation functions etc, in terms of the Askey- Wilson 
polynomials was achieved without any reference to the AW algebra. The ultimate relation of 
the exact solution in the stationary state to the AW polynomials was possible due to the AW 
boundary hidden symmetry of the ASEP with most general boundary conditions. The relation 
of the AW algebra to the if- matrix, determined by (63) and (64), and satisfying the reflection 
equation, puts a solution beyond the stationary state into perspective. 



7 Interpretation of the ASEP boundary operators 

As known the open ASEP is related to the integrable spin 1/2 XXZ quantum spin chain 
through the similarity transformation F = —qU~^HxxzU^ (see [13] for details). Hxxz is the 
Hamiltonian of the Uq{su{2)) invariant quantum spin chain (57) with anisotropy A^ and with 
added non diagonal boundary terms Bi and B^. 

Hxxz = H^xz + Bi + Bl (93) 

The transition rates of the ASEP are related to the boundary terms in the following way (/i is 
a free parameter, irrelevant for the spectrum) 

Bl = — (a + 7 + (a - -f)al - 2a/i(Tf - 27/i^ Vi+) (94) 

2g 



Br 



2g 



It has been shown by Sandow and Schuetz [56] that the bulk driven diffusive system with 
reflecting boundaries can be mapped to the spin 1/2 ?7q(sM(2))-invariant quantum spin chain. 
The Uq[su{2)) generators satisfying eqs.(74) and (75) act on the tensor product representation 
space {V^)®^ as 



A± = ^^q'^ ® ... ® q'^ ® af ® q^~^ ® ... ® g~ 



q^^ = q^ 2 ^q^2 ®...®q^2 (95) 

3. ZA A- ^ZA '"3 



where 173,0"^ are the Pauli matrices and the index i means that the matrix is associated with 
the ith site of the chain (ith position in the tensor product). The representation is completely 
reducible, the product of L spin 1/2 representations decomposes into a direct sum of spin j 
irreducible representations with maximal highest weight j = L/2 decreasing by 1 to j = 
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or j = 1/2 for even L or odd L. Within the matrix product approach the bulk process with 
reflecting boundary conditions is described by a quadratic algebra 

DiDo - qDoDi = (96) 

which defines a two-dimensional noncommutative plane with the SUq{2) action as its symmetry. 
The operators associated with the bulk ASEP form the two- dimensional comodule of SUq{2). 
As a consequence of eq.(96), for generic q, a spin j representation of Ug{su{2)) can be realized 
in the space of the g-symmetrized product of L = 2j two dimensional representations D^, 
/i = 0, 1, with basis Dq^'^D'i, k = 0,1,.. .,L. The stationary probability distribution, i.e. the 
ground state of the Ug{su{2)) invariant Hamiltonian H^-^^^ corresponds to the g-symmetrizer of 
the Young diagram with one row and L boxes [57]. The presence of the boundary processes (i.e. 
the nondiagonal boundary terms in the Hamiltonian) reduces the Uq{su{2)) bulk invariance and 
amounts to the appearance of linear terms in the quadratic algebra. The boundary conditions 
define the boundary operators which carry a residual symmetry of the process. It is expressed 
in the fact that the boundary operators are constructed in terms of the Uq{su{2)) generators, as 
seen from the explicit formulae (76). With A±^N being the generators of a finite dimensional 
Uq{su{2)) representation, it can be verified from eq.(76) that aDo—'jDi commutes with H{q)'^'^^ 
and f3Di — 6Dq commutes with H{—q~^)'^'^^, where according to [30] 



ffQGr^_^-i^ = -UH'^^^{q)U-^ (97) 

and 

U = exp\i-y] ma'i 1 (9^ 




Thus the boundary operators constructed as the linear covariant objects of the bulk Uq{su{2)) 
symmetry acquire a very important physical meaning - they can be interpreted as the two 
nonlocal conserved charges of the open ASEP. Such non local boundary symmetry charges 
were originally obtained for the sine Gordon model [9] and generalized to affine Toda field 
theories [TU] , and derived from spin chain point of view as commuting with the transfer matrix 
for a special choice of the boundary conditions [TT] . In particular, the left boundary operator 
aDo — 'yDi in the finite dimensional representation (76) is analogous to the one boundary 
Temper ley-Lieb algebra centralizer in the "nondiagonal" spin 1/2 representation [T2] . 



8 Discussion and conclusion 

In this paper we have considered the homomorphism of the Askey- Wilson algebra with two 
generators into the quantum affine Uq{sl{2)) algebra. This homomorphism defines the Askey- 
Wilson algebra as a coideal subalgebra of Uq{sl{2)). We have constructed an AW operator 
valued K matrix which is a solution of the boundary Yang-Baxter equation (refiection equa- 
tion). We consider the relation of an AW algebra to a solution of the refiection equation to be 
important for the exact solvability of a physical system with the quantum Uq{sl{2)) invariance 
in the bulk and hidden boundary Askey- Wilson algebra symmetry. 



As an example of a physical system with boundary AW algebra we consider a model of 
nonequilibrium physics, the open asymmetric exclusion process with general boundary condi- 
tions. This model is equivalent to the integrable spin 1/2 XX Z chain with nondiagonal bound- 
ary terms whose bulk invariance (infinite spin chain) is Uq{su{2)). The presence of boundaries 
breaks the bulk quantum affine symmetry of the equivalent quantum spin chain, however a 
remnant of the bulk symmetry survives and it is expressed in the possibility to construct the 
ASEP boundary operators in terms of the Uq{su{2)) generators. Thus the boundary operators 
of the open asymmetric exclusion process generate an Askey- Wilson algebra, which is the hid- 
den boundary symmetry of the process. The exact solution of the ASEP with most general 
boundary conditions (four boundary probability rates) in the stationary state was obtained [SS] 
in terms of the AW polynomials without reference to the AW algebra. It was emphasized that 
the solution was ultimately related with the AW polynomials. Such an ultimate relationship 
is natural from the point of view of the boundary AW algebra. The existence of the reflection 
matrix K{z,p) (and its dual K^{z~^,p*)) constructed in terms of the AW algebra generators 
and satisfying the boundary Yang-Baxter equation is, in our opinion, the deep algebraic prop- 
erty of the open asymmetric exclusion process that may allow for extending its exact solvability 
beyond the stationary state. 

It is important to emphasize the representation dependence of the Askey- Wilson algebra 
(as well as of the MPA bulk quadratic algebra (68)). Constructed as a coideal subalgebra, it 
has the property that the structure constants p, p*,uj,ri,ri* carry the information of the cor- 
responding quantum algebra Uq{su{2)). The boundary Askey- Wilson algebra whose structure 
constants depend on the finite dimensional Uq{su{2)) representations is the ASEP hidden sym- 
metry and this may have an important consequence in relation to Bethe Ansatz integrability. 
The Bethe solution of the open ASEP [SH] was achieved through the mapping to the Uq{su{2)) 
integrable XX Z quantum spin chain with most general non diagonal boundary terms, provided 
a particular constraint on the model parameters was satisfied. Quite surprisingly the constraint 
coincides with the condition for a finite-dimensional representation of the Askey- Wilson bound- 
ary algebra. The suitably chosen representation dependent boundary algebra may turn to be 
the key in relation to Bethe Ansatz integrability. For the ASEP the reduction of the bulk 
invariance gives rise to the boundary symmetry which remains as the linear covariance algebra 
of the bulk Uq{su{2)) symmetry and one can further employ Bethe ansatz to obtain exact re- 
sults for the approach to stationarity at large times and to completely determine the spectrum 
of the transfer matrix. As commented in [SH] the way one can satisfy the condition for the 
Bethe Ansatz solution of the ASEP implies additional symmetries. In our opinion, the linear 
covariance Askey- Wilson algebra of the bulk Uq{su{2)), whose generators are interpreted as 
the two nonlocal conserved charges of the ASEP, is the hidden symmetry behind Bethe Ansatz 
solvability. 

It is worth mentioning that the relation of the ASEP (or the equivalent quantum spin 
chain) boundary algebra to Bethe Ansatz integrability is promising form the point of view of 
Bethe Ansatz perspective in string theory. One is interested in closed strings with periodic 
boundary conditions. However, it is simpler to find the scattering matrix on the infinite line 
using asymptotic states and bootstrap. Then the spectrum is determined by asymptotic Bethe 
equations (HDIEI] and they are approximate for a system of a finite size. The study of the Askey- 
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Wilson algebra of a system on a ring with periodic boundary conditions which is interesting on 
its own, might be also useful for application to strings of finite length. 

We have obtained an Askey- Wilson algebra as a coideal subalgebra of the quantum affine 
Uq{sl{2) and implemented it to find a solution of the reflection equation. We have related this 
consideration to a model of nonequilibrium physics where the boundary operators generate a 
tridiagonal Askey- Wilson algebra, which is the linear covariance algebra of the bulk Uq{su{2)) 
symmetry. It is the hidden symmetry that allows for the exact solvability in the stationary state 
and provides the framework for employing Bethe ansatz to determine the dynamical properties 
of the open process. 
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